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Abstract 

The stability properties of line solitary wave solutions of the (2+l)-dimensional Boussinesq equa- 
tion with respect to transverse perturbations and their consequences are considered. A geometric 
condition arising from a multi-symplectic formulation of this equation gives an explicit relation 
between the parameters for transverse instability when the transverse wavenumber is small. The 
Evans function is then computed explicitly giving the eigenvalues for transverse instability for 
all transverse wavenumbers. To determine the nonlinear and long time implications of transverse 
instability numerical simulations are performed using pseudospectral discretization. The numerics 
confirm the analytic results, and in all cases studied, transverse instability leads to collapse. 

PACS numbers: 05.45.Yv, 47.35. +i 
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I. INTRODUCTION 



One of the fundamental ways that a solitary wave traveling in one space dimension 
generates a two space dimensional pattern is through transverse instability A transverse 
instability of a line solitary wave is associated with a class of perturbations traveling in a 
direction transverse to the basic direction of propagation. In addition to establishing the 
existence of transverse instability, a major question is what implications this instability 
means for the long-term behaviour of the system: does it settle into a new two-space- 
dimensional pattern, or collapse ? In this paper we study this sequence of questions for the 
canonical Boussinesq equation in two space dimensions 

u t t = (f(u) + eu xx ) xx + auy y , (1) 

where e = ±1 and a = ±1. In general, f(u) can be any smooth function, but the canonical 
form of the Boussinesq equation has the form 

f(u) = D{u 2 - u) with D = ±1 . 



When D — — 1, e — 1 and a = 1 this equation was derived by Johnson ||T8| to describe 
the propagation of gravity waves on the surface of water, in particular the head-on collision 
of oblique waves, and it was derived by Breizman and Malkin M in the context of Langmuir 
waves. 

In the absence of the transverse variation (i.e, u y = 0) and for e = —1, D = — 1 this 
equation reduces to the so-called "good" Boussinesq equation, which is well-posed, and for 
which sech 2 -solutions exist for any c with |c| < 1. These waves are stable when ~ < |c| < 1 
||. For the case |c| < \ it was shown by computer-assisted simulation of the leading term 
in the Taylor expansion of the Evans function that there is an unstable eigenvalue 0. This 
result was generalized to include solitary waves with nonzero tails, and rigorously proved 
using the symplectic Evans matrix in [[13] . 



Transverse instability of solitary waves has been widely studied since the seminal work 
of Zakharov [|25| on the nonlinear Schrodinger equation and the work of Kadomtsev & 
Petviashvili [19[ on transverse instability of the Korteweg-de Vries soliton. Since then, 
transverse instability of solitary waves has been investigated for a wide range of models; 



examples include the nonlinear Schrodinger (NLS) equation and related equations pTl [22 
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24| , Kadomtsev-Petviashvili equation 0, [5], [17|, |2(J, the Zakharov-Kuznetsov equation 



10, 22], and water waves [11]. A review of transverse instability for NLS and other related 



models can be found in Kivshar & Pelinovsky p0 |. 

In this paper, we will first use a geometric condition as derived in [jTIJ to get an explicit cri- 
terion for small transverse wavenumber instability. For this we use the multi-symplectic for- 
mulation of (|I]) in an essential way. To get detailed information for all transverse wavenum- 
bers we compute explicitly the Evans function for the (2+l)-dimensional Boussinesq model 
linearized about a larger family of line solitary waves (allowing the state at infinity to be 
nonzero). Plots of the dependence of the growth rate on the transverse wavenumber are 
presented. 

The post-instability behaviour of the nonlinear problem is studied using direct numerical 
simulation. The numerical evidence confirms the analytic results and suggests that the post- 
instability in the nonlinear system leads to collapse in all cases. A multi-symplectic pseu- 
dospectral discretization [15] is used as a basis for the numerical simulations. The numerical 
scheme is applied to the full two-dimensional PDE and we observe transverse modulation 
and further development of the longitudinal and transverse instabilities, resulting in the 
collapse of the initial line solitary waves. In the parameter region where the analytic crite- 
rion indicates that the solitary wave state is longitudinally stable but transversely unstable, 
simulations support the analytic results and provide insight into the long-term development 
of this instability. 

II. MULTI-SYMPLECTIFYING THE EQUATIONS 

The Boussinesq system has a range of geometric structures. Firstly, we record the La- 
grangian and Hamiltonian structures. Let u = <fr xx , then the system is Lagrangian with 

L = J [-\€t + F{<Px X ) + \e<Pl xx + \a4> 2 xy ]dxdydt, 

where F(-) is any function satisfying F'(-) = /(•). 

The Boussinesq equation can be represented as a Hamiltonian system in a number of 
ways (e.g. p3|). For example, let 

H = J [F(u) - \eu 2 x + \§l + \ow\ + 7(u - w x )] dxdy , 
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where 7 is a Lagrange multiplier associated with the constraint u = w x . With Hamiltonian 
variables ($, u, w,j) the governing equations take the form 



-u t 



xx 



8H 

8u 



= g 

= ^ 

07 



f(u) + eu^ + 7 

CTW„ 
M — W 7 



(2) 



7x " ) 



However, the most interesting form of (HI) for the present purposes is the multi-symplectic 



formulation which can be represented in the canonical form [14 

MZ t + KZ X + LZ y = VS{Z) z em 



where 





MZ t 














M = 




Pi 




P2 










( ° ' 




( 




( 


K = 








-1 ( 




- 




V ' 



(3) 



1 o\ 

-1 









) 1 o\ 

) 10 
) 1 
) 
1 
} -1 / 



with u(x,y,t) = qx(x,y,t) , 
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\ 00000/ 



S(Z) = -F( qi ) - ip? + \ V \ - a - V \. 



Using qi — u it is straightforward to show that this system is a reformulation of ([I]) . 



III. GEOMETRIC CRITERION FOR TRANSVERSE INSTABILITY 



An advantage of the multi-symplectic formulation is that there is a geometric condition 
which is easy to verify for transverse instability of line solitary waves [1(J . 
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Consider the well-known basic family of solitary waves of (|l|) of the form 

Z(x, y, t) = Z(8; c, /), 6 = x — ct + ly + 8 , 
obtained by taking the first component to be a sech 2 wave, 

u(6; c, I) = (e^ Z(6; c, I)) = A(c, I) sech 2 (£(c, 1)6) , 

with 



(4) 



(5) 



B{c, I) = \y/e{D + c 2 - aP) , A(cJ) = Q-B 2 . 

Existence of the solitary wave clearly requires e(D + c 2 — al 2 ) > 0. The other components 
of Z are easily obtained from (|5|) and the multi-symplectic equations (|3|). 

For the linear stability analysis, let Z(x, y, t) = Z(8; c, I) + $l[U(6; A, k)e xt+lky ], substitute 
this into ([3]) and linearize. Then, if the resulting linear equation has square-integrable 
solutions U(6; A, k) with 3ft(A) > and k G K., we call the basic solitary wave state Z(6; c, I) 
transversely unstable. Assuming that Zg is the only square integrable element in the kernel of 
the linearization operator C = D 2 S(Z) — [K — cM + ZL] ^, we have the following geometric 
condition of transverse instability for small A and k. Suppose 



A 



> 0, where 



A=-lfZ (MZ e ,Z)d9, 
B=^rj L Z e ,Z)d6, 



(6) 



Then the basic solitary wave Z(6; c, /) of (0) is linearly transverse unstable [TU, |TT|| . 
Using the above definitions of the multi-symplectic matrices M and L, we obtain 

1 f°° ( d d \ f°° 

A = '2 L [ qi de q2 ~ q2 de qi ) de = - c L qldd = ~ cK ' 

1 f°° ( d d \ f°° 

B= 2j^ [ qi de q3 ~ q3 de qi ) de = al L qldd = alK > 



(7) 
(8) 



where 

4 A 2 _6e 
3~B ~ ~D~ 2 

Substitution of (|7|) and (§) in (^) yields: 



K 



u 



2 de =-— = -^( a l 2 -c 2 -D) 



I 2 — c 2 — D 



a 



sgn 



sgn A = sgn 



a [ K + c—K + l-K 









-~ c A (a c + l -A^j 


= sgn 


a \A C + l -Aij 









sgn [—a [al 



2 c 2 



D) (4a/ 2 - Ac 2 - D)] . (9) 
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Since the condition for transverse instability requires A > 0, we have the following result: 

Suppose 

so- (4a/ 2 - 4c 2 - D) > 0, (10) 

then the basic solitary wave Z(9; c, I) is linearly transversely unstable. 

The multi-symplectic formulation also provides an expression for the linear growth rate 
of the instability A as a function of the transverse wavenumber k for long-wave perturbations 

13: 



This provides the growth rate for k small. In the next section, the Evans function will be 
constructed in order to determine the growth rate for all transverse wavenumbers k. 

In the remainder of this section, we apply the condition (|T0|) for various parameter values. 

For the "good" Boussinesq equation from || with e = —1 and D = — 1 the existence and 
transverse instability requirements are 

at 2 - c 2 + 1 > and -a (^al 2 - c 2 + ^ > 0. (12) 

respectively. Combining these conditions leads to the following system of inequalities for c 
and I when a > 

X - + al 2 < C 2 < 1 + al 2 , (13) 



and for a < 



c 2 < \ + ol 2 (14) 



These inequalities define the regions in (c, I) parameter plane, where the basic solitary wave 
exits and is linearly transversely unstable, and these regions are presented in Figure 1. 
One can do a similar analysis for Johnson's equation [|l8j, where a = 1, e = 1 and D = —1. 



The existence requirement is I 2 < c 2 — 1 and the instability condition is I 2 > c 2 — \. This 
result is inconclusive for two reasons. First, the two regions do not overlap so the geometric 
condition does not predict instability for any parameter values. Secondly, when e — +1 
the equation is ill-posed as an evolution equation (this can be seen at the linear level where 
the dispersion relation predicts instability as the wavenumber goes to infinity), and so the 
question of long time stability is irrelevant. 
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FIG. 1: Theoretical boundaries of transverse instability, a) Case ( |l~3|) with a = 1. The waves are 
unstable for the parameters lying within the shaded regions, b) Case (|14|) with a = —1. The waves 
are unstable for the parameters within the circle. 

IV. EVANS FUNCTION ANALYSIS OF THE TRANSVERSE INSTABILITY 

In this section we use the Evans function formalism in order to analyze the linear trans- 
verse stability problem for the Boussinesq model ([I]) for all values of the transverse wavenum- 
ber. We restrict attention to the parameter values of most interest: s — — 1 and D = — 1 
associated with the "good" Boussinesq, although we put no restriction on a (but keeping in 
mind that a = +1 is the most interesting case). 

However, the class of solitary waves will be enlarged. Namely, we include solitary waves 
bi-asymptotic to a nontrivial state at infinity, specifically, 

U(0) = U OQ + Q5 2 sech 2 (59) , 9 = x - ct + ly, (15) 

where 

6= - JVTT^-c 2 + aP and = |(1 - c 2 + al 2 ) - 25 2 = 2 ° (16) 

2 v 2 1 + y/1 + 4a 

The value of the parameter a is constrained only by existence of the square root: 1 + 4a > 
(c 2 - al 2 ) 2 . 

Here we will not use any geometric structure (although it might be interesting to look 
more closely in this direction) and so work directly with (|l|). Let 

u(x, y, t) = U(9) + &(u(0) exp [iky + At]). (17) 
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By substituting this expression in ([I]) and linearising, one obtains the following equation for 
the complex function u(9) 

ueeee + 2(Uu) 00 - (1 - c 2 + al 2 )u ee - 2(cA + iakl)u e + (A 2 + ak 2 )u = 0. (18) 

After the change of variable x = 56, substitution of the explicit expression for U from (|H>|). 
and dropping the tildes, equation ( P~8| ) reduces to 

Uxxxx - 4 [(1 - 3 sech 2 x)u] xx - ju x + (3u = 0, (19) 

where 

2(c\ + iakl) , „ \ 2 + crk 2 
7 = ^^ Und(3 = ^^. (20) 

To obtain explicit solutions of this equation, we note that by taking u = <p xx and v = 
1 — 3 sech 2 x in (|T9"D, and integrating twice the equation simplifies to 

<f>xxxx ~ 4v(f) xx - 70a, + p<j) = 0. (21) 

Solutions of this equation can be readily found in a manner similar to that in (see also 
0). First we note that in the limit x — > ±oo, equation fl2~lD reduces to 

<t>xxxx ~ ^xx ~ J(frx + (3(f) = 0. (22) 

Substituting now = e^cj), one can see that \i satisfies the quartic equation 

/i 4 - 4/i 2 - 7/i + /? = . (23) 

Quartics of this form have been analyzed in ]nj (see equation (10.9) there), and when 



3?(/3) > there are two roots with positive real part and two roots with negative real part. 
Therefore, the space of solutions decaying as x — > +oo is two-dimensional, as is the the 
space of solutions decaying as x — > — oo. 

If the four roots fij, j = 1,..,4 of the equation ( p3|) are distinct, the corresponding 
solutions of (pll) are given by 

^.( x ) = e ^ x hj(x), (24) 

with 

/^•(x) = (4^J + 8/Xj - 7) - 12/i 2 tanhx. (25) 
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The case of multiple roots can be handled similarly ||. Solutions of the original equation 
(|T9|) are found by substituting u(x) = 4>(x) xx , and the other components of the vector v(x) 
can be obtained by the differentiating the expression for u(x). 

Localised solutions of the linearised problem exist if one can match the solutions decaying 
as x — > oo with the solutions decaying as x — > — oo. This can be determined by finding the 
zeros of the so-called Evans function which correspond to the eigenvalues of the linearised 
problem. To define the Evans function, we write the equation (fl9| ) as a first-order system 

/ u \ ( 1 \ 



v x = A(x)v, V 



u x 



A[x) 



1 
10 

1 

\-P + Av xx 7 + 8v x Av J 



(26) 



with v — 1 — 3 sech x. 

Since the trace of the matrix A(x) vanishes, the Evans function can be defined as 
E(X, k) = Vx(x) A v 2 (a;) A v 3 (a;) A v 4 (a;) 0. An alternative expression for the Evans function 



can be derived by using the adjoint system as shown in ||12|| . The adjoint system of (|26|) has 
the form: 

f Wl \ 



w x = -A(x)*w, w 



w 2 
\w 4 J 



(27) 



where A(x)* denotes the Hermitian conjugate of A (A(x)* = A(x) ). The equation for w 4 
turns out to be 

(wi)xxxx ~ ^v{wi) xx + 7(w 4 )x + P(w 4 ) = 0. (28) 

This equation is equivalent to ( pTf ) up to the change of variables: x — > —x, 7 — > 7, j3 — * ft, 
and therefore its solutions can be obtained from (|24]) by changing x for — x and conjugating 
them: 

( W4 ) . = e -»3 x h£^j, (29) 
with hj(x) defined in (f25l). Other components of the vector w(x) can be obtained from 



Let Hi and [i 2 be the two roots of the equation (|23|) with negative real part, and let Vj(x) 
and Wj(x), j = 1,2 be the corresponding solution vectors of the linearised (respectively, 
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FIG. 2: a) The Evans function E(X) = E(X,0) versus A for c = 0.25, c = 0.35 and c = 0.75 
respectively, b) Growth rate versus transverse wavenumber for the values of velocity c = 0.6, 
c = 0.75 and c = 0.9 respectively. 

adjoint) system. Since the matrix A(x) in fl26|) is traceless, we can define the Evans function 
for the system fl26|) as follows [fP^] : 



E(X,k) 



(wi(0),vi(0)) (wi(0),v 2 (0)) 
(w 2 (0),v a (0)> (w 2 (0),v 2 (0)> 



(30) 



where (■, •) denotes the complex inner product in C 4 . To obtain a unique definition of the 
Evans function, the scaling \im x ^ 00 e~ 2 ^ jX ('Wj(—x),\j(x)) = 1 is used. This normalises the 
eigenvectors and the adjoint eigenvectors of A°° = lim^ioo A(x). 

After some lengthy algebra and introducing the scaling, which enforces the asymptotic 
limit E(X, k) — > 1 as A — > oo, the final expression for the Evans function can be obtained, 
which we do not present here since it is lengthy (the expression for the Evans function as 
well as the calculations of the instability growth rate can be downloaded as a MAPLE-file 
from the website 0). 

Zeros of the Evans function E(X, k) correspond to the bounded solutions of the linearised 
stability problem with the wavenumber k and the growth rate 3?(A). The leading order 
terms (in k and A) in the Evans function are in complete agreement with the results of the 
geometric condition of §3. Note that, since the construction here is based on a basic solitary 
wave with a nontrivial state at infinity, it is suggestive that the geometric condition [T^ 
extends to such waves. 
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We illustrate the dependence of the Evans function on the wavespeed and transverse 
wavenumber in Figure 2. In the left graph, the transverse wavenumber is set to zero, to 
compare with known results on longitudinal instability. The graph is in complete agreement 
with known results (e.g. [0, that the solitary wave is stable for | < c < 1 and unstable for 
< c < |. In the right-hand graph in Figure 2 we present the plot of the growth rate 3?(A) as 
a function of the transverse wavenumber. Note that waves of the good Boussinesq which are 
longitudinally stable are transverse unstable. Note also that there is a cut-off wavenumber, 
similar to other cases of transverse instability, such as in the Zakharov-Kuznetsov equation 



V. POST-INSTABILITY SIMULATIONS 

In this section we perform a simulation of the PDE (H) using the multi-symplectic spectral 



1 2N 



discretization proposed in [[jjj and applied there to Zakharov-Kuznetsov and shallow-water 
equations. 

The (2+l)-dimensional Boussinesq equation is considered with e = D = — 1 on a finite 
domain (x,y) = [0, L\ x [0, L] with L > some constant, and periodic boundary conditions 
on both spatial variables. We choose a spatial mesh-size as Ax = Ay = Am = L/2N and 
introduce the discrete two-dimensional Fourier transform defined as 

2N 

2X 



where 

e k - 1 — 

and Uij ~ u(mij),mij — (i — l)Ax + (j — l)Ay (cf. ]TB[). Fourier spectral discretization 
of the (2+l)-dimensional Boussinesq equation yields 

duU kl = e 2 k [e6 2 k U kl + VjwF(U)] + «T0fE4i, (31) 

where 9k are the entries of the diagonal matrix defined by the relations 

9 k = 9 k , for k= 1,...,N, 

9n+i = 0, and 

h = -02N-k+2, for k = N + 2, 2N 
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FIG. 3: a) Development of the longitudinal instability and collapse at time t = 12 for c = \. b) 
Propagation of a stable solitary wave for c = § . 



FIG. 4: Energy evolution. Dashed line represents the initial energy level, and the solid line shows 



r 



the time evolution of energy, a) Unstable case c = \ . b) Stable case c 



which follow from the periodicity of the discrete Fourier transform [jnj, and F(XJ) denotes 
the Fourier transform of the anti-derivative of the function f{u) in (|l|). The same result would 
be obtained if one applied the spectral discretization to the multi-symplectic formulation 



as it was done for the Zakharov-Kuznetsov equation in [15 



For the second-order time derivative we used the central difference approximation (time 
step was chosen to be At = 0.01 in all the simulations): 



duU, 



U^ +1 - 2U% + uir 1 



At 2 



(32) 



One should note that the only valid test of this scheme can be done for the "good" Boussinesq 
equation with a > 0. For o < in the case of the "good" Boussinesq equation, an initial 
profile independent of x would result in a solution which could grow "faster than exponential" 
because for large transverse wavenumbers, growth rate of the initial data has no upper bound 
(ill-posedness). 

To test the algorithm, we first used it to confirm the results for the dynamics of the 
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FIG. 5: Solitary wave for a = 1 and c = \. a) Initial profile, b) Development of the transverse 
modulation (time t = 11.25 ). 




FIG. 6: The same as in Fig. 5. a) Wave collapse ( time t = 12). b) Energy evolution. Dashed line 
represents the initial energy level, and the solid line shows the time evolution of energy. 



one-dimensional solitary waves. The initial profile was taken to be of the form 



u\x) 



°- (1 - c 2 ) sech 2 



i(l-C) 



L 

X ~2 



(33) 



where £(x) is a small random perturbation. The results are presented in Figures 3 and 
4. For c = \ the solitary wave solution is linearly unstable as reported in || 13| , and the 
development of this linear instability is shown in Figure 3a). In the case c = | the numerical 
results confirm the stability of the solitary wave (see Figure 3b) ). The simulations were 
run on an interval of the length L = 256 with 2N = 512. As a numerical check, the total 
energy was monitored, and it was found to be well behaved till near the collapse when the 
significant errors occur, as illustrated in Figure 4. 

For the two-dimensional simulations we took an initial profile in the form of the line 
solitary wave uniform in y 

i 1 — ? I scrlr ^ 



u(x, y, 0) = — (l — c 2 ) sech 2 



(34) 



where £(x, y) is a small random perturbation (in this case 1 = 0). The length of the square 



13 




y 



FIG. 7: Solitary wave for a = 1 and c = |. a) Initial profile, b) Development of the transverse 
modulation (time i = 184.5). 

a) t>) „. , , , , . 




FIG. 8: The same as in Fig. 7. a) Wave collapse (time t = 246). b) Energy evolution. Dashed line 
represents the initial energy level, and the solid line shows the time evolution of energy. 

box was chosen to be L = 64 with the number of Fourier modes 2N = 128. In the case 
c — j the solitary wave flP|) is linearly unstable in longitudinal direction as is known from 
the stability analysis of the ID equation. In Figure 5b) we can see this instability developing 
in a similar way as in the ID case. One can also note in this Figure the development of 
the stable transverse modulation. Wave collapse in this case is shown in Figure 6a), with 
the plot of energy as a function of time in Figure 6b). When c = |, the solitary wave is 
longitudinally stable but transversely unstable, and the development of this instability is 
presented in Figures 7 and 8. We note that at the initial stage of the evolution there is 
a transverse modulation developing while the amplitude of the wave is gradually growing 
Figure 7b), then the instability prevails leading finally to the collapse of the wave Figure 8a). 
Note that this collapse is clearly a two-dimensional effect, since it does not happen uniformly 
in the ^/-direction. The energy proves to be conserved rather well during the simulations (see 
Figure 8b) ), although the energy deviates substantially as the wave approaches the stage 
of collapse. 
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VI. CONCLUDING REMARKS 



We have considered the transverse instability of line solitary wave solutions of the (2+1)- 
dimensional Boussinesq equation. Using the multi-symplectic formulation of the system, we 
derived a geometric condition for this instability for small transverse wavenumbers. With an 
Evans function approach, the linearised stability equation was analyzed, and this allowed to 
obtain the dependence of the instability growth rate for all transverse wavenumbers. Numer- 
ical simulations support the analytical results about transverse and longitudinal instabilities 
and demonstrate the development of those instabilities and subsequent wave collapse. 

We conclude with an open problem. While analytic theories for collapse of solitary waves 
for the Boussinesq equation in one space dimension exist ||23|| , it is an interesting open 
problem to develop an analytical technique for predicting collapse for the case of two space 
dimensions, e.g. a generalization of the virial theorem or the result of ]23| for example, and 
moreover, to determine if transverse instability for ([J) leads to collapse for all parameter 
values. 
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